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Techniques used to explain the nature of the optical effects of clouds on the light produced 
by lightning include a Monte Carlo simulation (10), an equivalent medium approach (6), and 
methods based on Boltzmann transport theory (4, 7). A cuboidal cloud has been considered 
using transform methods (7) and a diffusion approximation (4). 

Many simplifying assumptions have been used by authors to make this problem tractable. 
In this report, the cloud will have a spherical shape and its interior will consist of a uniform 
distribution of identical spherical water droplets. The source will be modeled as a Hertz dipole, 
electric or magnetic, inside or outside the cloud. An impulsive source is used. Superposition 
may be employed to obtain a sinusoid within an envelope which describes a lightning event. The 
problem is investigated by transforming to the frequency domain, obtaining Green’s functions, 
and then using the Cagniard-DeHoop method (3) to symbolically recover the time domain 
solution. 

The Hertz dipole potential may be derived from a wave equation whose source is propor- 
tional to the dipole moment per unit volume, 

P = Np, p = Qs, [1] 


where N is the number of dipoles per unit volume and p is the dipole moment of a single 
dipole constructed from equal and opposite charges of magnitude Q separated by a distance 
s oriented along s. From Maxwell’s equations, it is possible to show that the electric Hertz 
vector is the solution to 


v 2 n - [ic 


d*n 

dt* 



[ 2 ] 


The resulting electric and magnetic fields may be written in terms of II by 

H = £ v A f, E = w.n-^. 


[3] 


Similarly, for a magnetic dipole, 

v 2 n m - = -m, E = /iVA^, h = vv-n ffl -^. [4] 

The symbols are defined as in (9). Polarizations may be induced in a dielectric by the fields. 
The dipole moments used above are considered to be controlled only by the source. The sources 
will be taken to be unit pulses (polarizations per unit time), 


P = £(r - r')S(t - t')p, M = <$(r — r ')6(t — f # )p m - 


[5] 


Superposition may be used to construct fields due to general temporal and spatial dependencies. 

In the following, script variables will denote quantities in the frequency domain. They are 
related to the time domain variables through 

CO oo 

E(r, t) = J db,€(r, «)«-*■*-*'>, H(r, t) = jdojH( r, w)e- <to,(t -*' ) . [6] 

-OO — OO 


The transforms of the Hertz potentials are given by 


,1X 


*=—hW(KR)P, f m =—h^(KR)p m , /£>(*)=—, K*=u a ne, 

47 re 47r zx 


[7] 
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where R= |r- r'|, r is an observation point and r' is a source point. The electric and 
magnetic fields due to these dipoles are given by 

? = ' v > v R+ (~ K ’ - ^+^)[(p • 

w = ~ (}K - j)iVRr,i,)hi:m<R) 

where Z m and 7? m are obtained from H and S by replacing by an( j ihu by i K 

I he electric and magnetic field vectors may be written in terms of Hansen’s functions (5) as 

£ ■ = ^7 E«n [ M iU( A '> r)Mi„ n (A', r') + Ni„(A, r)Ni mn (A, r')] • p, 

v - v v ™ ( 2 »+i)(-ir . , 191 

X = E X ■ «»-— (n+1) ■ 


n=i m =—n 


for an electric dipole source; the magnetic field is obtained from £ by replacing —■ by 
and interchanging ML mn (A', r') and NL m „(A\ r'). The superscript 1 denotes a Hansen func- 
tion containing regular (Bessel type) functions; 3 refers to outgoing (Hankel type) func- 
tions. Analogous forms for fields due to magnetic sources are obtained by replacing the 
constants to the left of the summations by an< * 7T respectively and interchanging 

ML mn (A’, r') and Nl mn (A', r') again. If r < r 1 , the corresponding vectors in [9] are inter- 
changed. 

If the medium inside the cloud is specified by K(n, c) and that outside the cloud by 
k(fi 0 , «o) and the cloud is taken to be a sphere of radius a, the fields may be written as 

+ £'mi &e = 7? s + H l , Gm = 7?^ + Hm, where the superscript on 
the Q s indicates whether the Green’s function corresponds to the electric or magnetic field 
and the subscript indicates the type of source. The superscripts i and s refer to the source (z) 
or the scattered field (s). Green’s functions for external excitation (1), which were obtained 
by enforcing continuity of the tangential components of the fields across the boundary, may 
be easily modified to work for sources inside a sphere. Explicit results for an electric dipole 
source are given by 

r y c{mL(A', r)[M 3 _ mn (A', r') + a n ML mn (A', r')] 

+ N^ n (A',r)[Ni mn (A’,r') +6„NL mn (A',r')]}, r < /, 

£ = ^^{[Mi n (A’,r) + a n M^ n (A^r')]ML mn (A', r') [10] 

+ [N3 n (Jir,r) + 6 n N|i lB (Jir,r')]NL mn (if ) rO}, r' < r < a, 

k Y^Qn [c n M^ n (fc,r)Mi mn (A^rO+d n NL(A:,r)Nl mn (A%rO], t> Qj 
when the source is inside the sphere (r' < a) and 

' ^0”KM^.(*-,r)Mi m „(t,r') + b„N^(Jf,r)Ni mn (t,r')], r < a, ' 

E «” { r) + (.MU*. r)]Mi mn (*, r') 

« = ^| + [N^ n (t,r) + » n N| M ,(*,r')]Ni„„(t,r')}, a < r <r\ . [ 11 ] 

E 0" { r) [ML mn ( k , r') + c„Mi „ (*, r')] 

. +N3„„(t,r)[Ni m „(t,r') + »„Ni m „(t,r')]}. r > r’ 
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when the dipole is outside the cloud. The associated magnetic fields are obtained by re plac - 
in 6 ^7 and i)H(r-a)} and {i+(^ - i )H(a - r)}, 0 = Jp , 

where H is the Heavyside function, respectively and interchanging those M’s and N’s which 
are functions of r but not those which are functions of r'. Vector fields are recovered by inner 
multiplication with the source polarization from the right (i.e., Q — Q • p). If the source is a 
magnetic dipole, and are obtained from Q\ and by replacing and ^ 

by — - , and respectively and interchanging only those M’s and N’s which 
are functions of r'. The constants are given by 

^(Xjn)' ~ »)'’ " P^(xjn)' ' 

, _ ±(xh n y - *?(xj n y [ 12 ] 

^(xjn)' - Pfirt n/’ " ^ l (F^)n)'-/?^(xjn) ,, 

X — ka, y = A a, jn — jn(x ), 3n — jn(f)> = ^(x)» bn hji ^(?)* 


The constants a n and b n are obtained from c n and d n respectively by multiplying by (3 
and interchanging h n with J n , x with jc, and j n with f) n ; constants c n and D n are obtained from 
a n and b n respectively by multiplying by — 1 and interchanging h n with J n and f)n with j n . The 
abovementioned interchanges apply only to the respective numerators. Allowing r ~ oo and 
using asymptotic forms for the Hansen’s functions gives explicit forms for the far-field scatter- 
ing amplitudes g. Obtaining g for external excitation of a water droplet allows calculation of 
the bulk parameters of the cloud using the methods described in (6). 

A dyadic version of Helmholtz’s surface integral representation as given in (11) for rp = 
£ e or £ h , 


^(r, r /) = J dS(r 5 ){ V 5 AG(r-r s )] • [n A ^(r 5 , r') j 
s 

- nAG(r-r 5 ) • [V 5 A $(r„ r')] j = |g, V>} , 
<5(r— r s ) = ^ ( ! +^r) h ’ h = R = r - r i5 


[13] 


where V 5 operates on the variables associated with r 5 on the surface 5 of the scatterer and T 
denotes the transpose, may be used with the asymptotic form of h to obtain 


G = -^{i-rrjh, ^(r, r') = hg (r, r'), g - ,fcrr *, r') }. [14] 

Here, r, r s , and r' refer to vectors to the observation point, surface point, and source point 
respectively. Using Sommerfeld’s integral representation (12) for h, we write 



in [13], interchange the order of integration, and recognize g from [14] to obtain 

1>{r, r')= J dQ c e tklc t [g (r c , r')]. 

C 


[15] 


[16] 
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This spectral representation is valid for r greater than the scatterer’s projection on r. 

Substituting the form of g obtained from [10, 11] together with asymptotic representations 
for Hansen’s functions in [16] or, equivalently, using Sommerfeld type complex contour integral 
representations for the M’s and N’s (5) which are functions of r gives 

^( r > r ') = r')^ J dQ c e iki ’ c ' r S— mn (r c ), [17] 

C 

where R is M or N, S is B or C, and tjk= K ; A involves the coefficients in [12]. Rotating the 
coordinate system so that f is parallel to z, we obtain f c • r = r cos 6 C . In this form, the time 
domain solution may be written in a symbolic form by inspection (3). We obtain, for ct > R, 


^(r, r', t) = 


?(!)• ' 


/ [Pff ( ft ) ] 

_0 

[« 2 -(f) 3 ] t 


[18] 


When ct < R, rp — o. The sum is to be taken over all relevant combinations and the forms 
containing derivatives are to be expanded in a series and operate from left to right. 
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